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We study the effects of gluon saturation on pion production in high energy proton-nucleus colli- 
sions using the color glass condensate model. At high p±, we show that the px-distribution of 
gluons behaves as ~ l/p± in accordance with both conventional perturbative QCD calculations and 
experiment. Fragmentation of gluons into pions leads to a rapidity dependent depletion of pions 
relative to the conventional perturbative QCD predictions. We argue that these clear and system- 
atic differences provide a signal for the onset of gluon saturation which is accessible in upcoming 
experiments. 

PACS numbers: f 2.38.Aw, f3.87.Fh 



Since the pioneering work of Ref. [Q , the phenomenon 
of gluon saturation, occurring at very high gluon density, 
has been under intense study. One hopes that saturation 
would naturally tame the rise of the gluon distribution 
function at small values of Bjorken x and circumvent the 
subsequent violation of unitarity. 

A promising attempt to address the problem of satu- 
ration involves replacing the partonic fields by coherent 
classical fields. The justification is that for large occupa- 
tion numbers the dynamics are essentially classical. This 
program was initiated in Ref. and developed further 
in e.g. Refs. ||(§. Within this approach one is able to 
show that an intrinsic scale, the saturation scale Qs,a(x), 
is generated. This scale depends on both Bjorken x and 
atomic number, A, and roughly delineates the dense and 
dilute regimes. Throughout this work we shall refer to 
this approach as the color glass condensate (CGC) frame- 
work. 

The parametric dependence of the saturation scale 
has been extensively studied. Phenomenological models 
based on the idea of parton saturation have been success- 
ful in describing both the inclusive and diffractivc HERA 
data for F 2 || . According to this model, the dependence 
of the saturation scale on Bjorken x is 



(x) = (lGeV) 2 (^) A A 1 / 3 , 



(1) 



where A = 0.25 ... 0.3 and x = 3 • 10~ 4 . For the nu- 
merical calculations, we will use A = 0.3. Since Q 2 S A can 
be interpreted as the transverse density of partons, it is 
enhanced by considering larger nuclei at any given x by 
a factor of A 1 / 3 . As a result, one may expect to observe 
saturation effects better in nuclear collisions at currently 
available collider energies. Indeed, the present data from 
BNL-RHIC seem to be compatible with saturation 
based phenomenology [lO|-[T^]. These results, however, 
are far from being conclusive on account of numerous 
model dependencies related to our poor knowledge of 



the details of the spacetime evolution of the produced 
strongly interacting matter. 

The onset of saturation has important qualitative and 
quantitative consequences for the p^-distribution of the 
produced gluons. This was already pointed out in Ref. 
|jj in the context of proton-proton collisions where an 
asymmetry in the saturation scales of the two protons ap- 
pears away from the mid-rapidity. As one hadron moves 
deeper into the saturation region, towards smaller x, its 
saturation scale increases, while the saturation scale of 
the second hadron decreases as is evident from Eq. ([!]). 
In proton-nucleus collisions, a similar situation is real- 
ized already at mid-rapidity due to the atomic number 
asymmetry. In particular, one expects three different re- 
gions in the transverse momentum distribution of the 
produced gluons on account of the asymmetries in the 
saturation momenta, Qs,A > Qs,p- There is the region 
V-L < Qs.p < Qs,a in which both the proton and the 
nucleus are saturated. The second possibility is that 
the nucleus is saturated but the proton is in the lin- 
ear regime (Qs,p < P± < Qs.a)- The final region is 
the usual perturbative domain in which both gluon dis- 
tributions are entirely governed by the linear evolution 
(Qs,p < Qs.a < P±)- Consequently, proton-nucleus col- 
lisions offer the best possibility of experimentally probing 
the systematics of gluon saturation. The problem of eval- 
uating gluon production in proton-nucleus collisions has 
been found to be analytically tractable above the smaller 
saturation scale and in Ref. [[uj this was qualitatively 
studied in the CGC framework. 

In this Letter, we undertake the quantitative compu- 
tation of the produced hadronic spectra based on the 
framework of Ref. and compare it to the conventional 
perturbative QCD (pQCD) calculation which does not 
include any saturation effects. The main emphasis is on 
a numerical calculation of the fragmentation of the gluon 
distributions into pions, and on the study of how satu- 
ration effects manifest themselves in the ^-distributions 
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of the measured hadrons. We shall see that although the 
spectra at the gluonic level are qualitatively very differ- 
ent, this difference is made less striking by the fragmen- 
tation of the gluons into pions. 

In pQCD, the inclusive cross section for gluon produc- 
tion is given as a convolution of the hard partonic scat- 
tering cross section with the gluon distribution functions. 
The hard cross section behaves as ~ but the fi- 

nal result receives additional contributions from the two 
gluon distributions and from the running of the strong 
coupling constant. As a result, the slope of the trans- 
verse momentum distribution actually falls much more 
rapidly with a scaling closer to ~ Any realistic 

model of particle production must reproduce this scaling 
at large transverse momenta. 

We consider a CGC model obtained by combining two 
separate limits of the results of Rcf . (m| . The gluon pro- 
duction cross section in the perturbative domain is given 
by 

^ g port 8N C (N 2 - 1) al 

= * ^Mx P ,P±)xa(xa,P ± ), (2) 

where 



irRi N 2 - 1 



\ { / dxg i (x,p 2 L ) + ^- [ dxq l (x,p 2 L ) ) . 



with i 



p±e 



TV 



A, p. Ri is the transverse radius, x^ p — 
s, Cf = {N 2 — 1)/ (2N C ) and the nucleus is cho- 
sen to have negative beam rapidity. Here qi and gt denote 
the valence quark and gluon distributions in the nucleus 
or in the proton. One needs to evaluate the color charge 
functions Xi only in the weak field regime, and therefore 
the use of DGLAP evolved parameterizations is appro- 
priate. 

In the saturation domain of the nucleus, the cross sec- 
tion is given by 



da s g at 
d 2 p±dy 



C{N* - 1) a 3 



2tt 2 



9 

PI 



Xp(x p ,p 2 j_)ttR 



A- 



(3) 



The constant C is chosen so that the sub-cross section is 
continuous at Qs,A- We have neglected the logarithms 
appearing in the results of Ref. |Q since one expects 
the overall behavior to be driven by the strong powers of 
p± and the constant under the logarithm is in any case 
difficult to fix. The overall normalization is unimportant 
as we will discuss later. 

From these equations, it is clear that retaining the 
DGLAP evolved parton distributions for the weak fields 
and letting the strong coupling run will affect the qual- 
itative results of Ref. [JlJJ which display the expected 
~ 1/Pj_ and ~ l/,Pj_ behavior in the perturbative domain 
and in the saturation domain of the nucleus, respectively. 



The inclusion of these effects is essential for making the 
proper contact with conventional result for high p± men- 
tioned above. 

In order to compute the produced hadronic spectra, we 
convolute the above inclusive gluon production cross sec- 
tion with the appropriate fragmentation functions. For 
the sake of simplicity, we do not attempt to account 
for the Cronin effect fl5|| , nuclear shadowing |l6| or nu- 
clear modifications of the fragmentation functions [fl7|| . 
So, xqa(x, Q 2 ) — Axg p {x,Q 2 ) in the above equations 
and similarly for quarks. The generalization to nucleus- 
nucleus collisions with A\ <C Ai is trivial through re- 
placements p — *■ Ai, A — ► A%. We note that the overall 
normalization of the results is difficult to fix. In Ref. 
]l8t it was found that in order to reproduce the data 
from hadronic collisions above a cutoff momentum of 
Po ~ 1 • • • 3 GeV using pQCD one needs a -^/s-dependent 
if-factor to account for the higher order corrections. The 
relative normalization between the pQCD and CGC re- 
sults, however, is almost fixed since both results should 
match at high p±. Therefore, we do not attempt to 
fine-tune the overall normalization, but rather consider 
a leading order pQCD calculation, including only glu- 
ons, without any iif-factors and point out the differences 
relative to the CGC calculation which will be normal- 
ized at high p± to the pQCD result. Our central results, 
the slopes of the p±- and y-distributions, are not sensi- 
tive to the normalization ambiguities. One should note, 
however, that the y/s growth of the gluon multiplicity 
in the CGC model is driven by the small x-growth of 
the gluon distributions, and is therefore of the order of 
dN/dy ~ y^ A - We will focus exclusively on pions since 
this distribution approximates well the distribution of all 
hadrons. 

To compute the fragmentation of gluons into pions, we 
assume that the pion is produced collinearly with its par- 
ent gluon, rjg = 7]^, and carries a fraction z of the parent 
gluon's energy. The cross section for the production of 
pions is 



dcrP A - 



d 2 q±dy 



A 7 drrP A ^9 x 



where 



q± 



P± = —J(rn_L,y), y g = sinh 



1 ' — — sinh y 

-1/2 



J(m±,y) 



and ni 2 ^ = q\ 



nri 1 ^ cosh y J 

; 2 r . The integration over z is limited 
by the maximal energy the gluon can carry and by the 
lower bound on p±. In pQCD, this latter scale is a fixed 
cutoff, while in the CGC calculation it is determined by 
the saturation scale of the proton, 
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am_i_ . / q 

— =- cosny < z < mm 1, 
\/s 



J(m±,y)J , (5) 



where a = 1, 2 for the CGC and pQCD calculations re- 
spectively. This difference arises from the underlying par- 
ton kinematics. In LO pQCD, one produces two minijets, 
back-to-back in the transverse plane, carrying at most 
E = i/s/2. In the CGC calculation, the scattering is a 
BFKL-type 2 — > 1 fusion, and the produced gluon can 
have energy up to E = y/s. We use the CTEQ5 param- 
eterization of parton distributions |ist | and KKP param- 
eterization of the fragmentation functions [^pf . 

In Fig. [|, we plot the pj_-distributions of the produced 
gluons at different rapidities. The topmost three curves 
are for y = and the lower two are for y = 3. The dashed 
curve is the CGC result including only gluons in the 
sources Xi ■ Note that at y — the slopes of the CGC and 
pQCD calculations almost match at high p±. Below the 
saturation scale of the nucleus, the two results strongly 
deviate. In a full computation of the gluon distribution, 
this sharp bend at Qs,a(v) will become smooth. The es- 
sential point is that there is a marked depletion of gluons 
for p ± < Qs,a{v)- 
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FIG. 1. Transverse momentum distributions of gluons pro- 
duced in a proton-gold collision computed from conventional 
pQCD (dotted) and from the CGC calculation (solid). The 
uppermost set of curves is for y = 0, while the lower ones cor- 
respond to y = 3. The dashed curve shows the CGC result 
with purely gluonic source. 

If this spectrum were measured, one should be able to 
determine whether the parton dynamics in the collision 
are better described by pQCD or by the CGC model. To 
obtain a measurable spectum one needs to fragment each 
of the distributions in Fig. [I] to hadrons. 

In Fig. |^, we show the pj_-distribution of the pions pro- 
duced from the two models. The fragmentation changes 
the qualitative behavior of the spectrum of the CGC 
model dramatically with the sharp bend disappearing 
completely on account of the redistribution of momentum 
from the gluons to the pions. We note that there are sig- 
nificantly fewer produced pions at lower transverse mo- 
mentum in the CGC calculation than in the conventional 
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FIG. 2. Transverse momentum distributions of pions pro- 
duced in a proton-gold collision computed by fragmenting 
gluons from conventional pQCD (dotted) and from CGC 
(solid) calculations. The topmost two curves are for y = 
and the other pairs are for y — 2 and y = 3, respectively. The 
dashed line shows the suggested fit of the pion distribution. 
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TABLE I. The fit parameters C(y) and b(y) in Eq. (g) at 
different rapidities. 

quantum pQCD one and that this deficit is strongly ra- 
pidity dependent. We have checked using a smooth fit 
interpolating between the two limiting forms of the gluon 
distribution that the qualitative, and to some extent even 
quantitative, results for the pion spectrum in the CGC 
model are not sensitive to the sharpness of the knee in 
the gluon distribution. 

To obtain an useful simple parameterization of the pion 
spectrum, we consider the form 



da 



C(y) 



dydpl {pl+ fi(y) 2 ) b (yy 



(6) 



where the emergence of a rapidity dependent scale, fi(y), 
is a natural consequence of the rapidity dependence of 
the saturation scale. The rapidity dependence of the ef- 
fective slope b(y) is expected since the slope of the gluon 
distribution already depends on rapidity. We find that 
M (jy) 2 = 0.50e Ay , where A = 0.3. The values of the fit pa- 
rameters C (y) and b(y) at various rapidities are collected 
in Table The inclusion of the quarks plays a role for 
the overall normalization and the slopes are left almost 
unchanged as can be seen already from Fig. [l] 

Finally, we consider the rapidity distributions at fixed 
p±. One cannot approach the beam rapidity of the nu- 
cleus in this framework, since at such large negative ra- 
pidities the saturation scales of the nucleus and the pro- 
ton become comparable and the underlying model breaks 
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FIG. 3. The rapidity densities at fixed p±. The dot- 
ted curves show the spectrum of gluons in the CGC model 
at p± = 1.4 GeV (upper) and p± — 3 GeV (lower). 
The solid curves give the pion spectra from CGC model 
at p x = 0.4,0.6,0.8,1.0,1.7 GeV (top to bottom), and the 
dashed ones show the pion spectra from pQCD at p± = 0.4 
(upper) and 1.0 GeV (lower). 

down. Near midrapidity, the CGC model gluon spectrum 
in perturbative domain is characterized by constant be- 
havior in rapidity, as is the case for pQCD spectrum. 
In the saturation domain, the CGC gluon spectrum has 
a nonzero slope also at midrapidity. The CGC model 
gluon spectra are shown in Fig. [| with dotted curves for 
p± = 3 and 1.4 GeV, and the corresponding pion spec- 
trum by solid lines for p± = 0.4, 0.6 0.8 1.0 and 1.7 GeV. 
The pion spectrum from pQCD, shown by dashed lines 
for 0.4 and 1.0 GeV has the flat behavior at midrapidity 
over the whole p± range. We observe that the small p± 
pions might reveal the characteristic form carried by the 
gluon spectrum in the saturation domain. This result 
is easy to understand, since the pion spectrum at scale 
q± effectively probes the gluon spectrum at some higher 
scale q±/z. For the perturbative distribution one finds 
that z~ 0.5. ..0.6. This is clearly seen by comparing pi- 
ons at p± = 1.7 GeV and gluons at p± — 3 GeV. As the 
slope of the gluon distribution decreases, z moves towards 
smaller values as well. This causes the gluons in satura- 
tion domain to fragment mainly to much smaller values 
of p± than the perturbative gluons would. The situation 
improves at LHC energies where the saturation domain 
extends already at midrapidity to ~ 2 GeV. At large ra- 
pidities it is difficult to distinguish between saturation 
and perturbative spectra, and presumably different dy- 
namics which are not contained in either of these models 
become dominant towards the fragmentation regions. 

We have investigated the quantitative behavior of the 
p±- and y-distributions of gluons and pions in collisions 
of nuclei with a large atomic number asymmetry. We 
have shown that within the color glass condensate frame- 
work one indeed reproduces the behavior of the conven- 
tional pQCD at highpj_ which significantly deviates from 
commonly cited expectation of 1 /p]_ in the literature on 
saturation. The fragmentation of gluons into pions was 



shown to significantly change the qualitative shape of 
the spectrum. However, a clear difference from the spec- 
trum obtained from the conventional pQCD calculation 
in the form of a rapidity dependent pion deficit can still 
be extracted by measuring the rapidity dependence of the 
p_i_-distributions. The rapidity dependence of the slopes 
were found to be different at small p±, and a change in 
the form of the rapidity distribution of pions as a function 
of p± could be used to determine the onset of saturation. 
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